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ABSTRACT 

We present a kinematic model of the neutral-hydrogen in the spiral galaxy of the lens system 
PKS1830-211, based on a Muhi Element Radio-Linked Interferometer (MERLIN) 1.4-GHz radio map 
and the integrated and redshifted 21-cm hydrogen absorption-line profile as measured with the West- 
erbork Synthesis Radio Telescope (WSRT). Degeneracies in the models do not allow a unique determi- 
nation of the kinematic center and forthcoming deeper Hubble Space Telescope observations with the 
Advanced Camera for Surveys (ACS) are required to break this degeneracy. Even so, we measure the 
inclination of the hydrogen disk: i = (17 — 32)°, indicating a close-to face-on spiral galaxy. The optical 
depth increases with radius over the extent of the Einstein ring, suggesting HI depletion towards the 
lens center. The latter could be due to star formation or conversion of HI in to molecular hydrogen 
because of a higher metalicity/dust content in the galaxy center. The neutral hydrogen optical depth 
gives iVni = 2 x 10^^ cm~^ at r = 5.0 h^^ kpc in the disk (Tg = 100 K), comparable to local spiral galax- 
ies. Our study shows that planned new radio telescopes (i.e. ALMA, LOFAR and SKA) - which will 
discover large numbers of similar lens systems - are powerful new tools to probe the internal structure, 
kinematics and evolution of spiral galaxies in detail to z > 1, as well as their neutral- hydrogen (and 
molecular) content, thereby complementing studies of HI emission from spirals at z < 1 and damped 
Ly-Q; systems to z 3> 1. 

Subject headings: gravitational lensing — galaxies: structure — ISM: kinematics and dynamics — line: 
profiles 



1. INTRODUCTION 

The gravitational lens PKS 1830-211 (e.g. Pramesh 
Rao & Subrahmanyan 1988; Subrahmanyan et al. 1990) 
is the brightest known radio lens in the sky. The source 
at z=2.51 (Lidman et al. 1999) is extended and lensed 
into an Einstein ring (e.g. Jauncey et al. 1991) by a 
high-redshift spiral galaxy at z = 0.89 (Wiklind & Combes 
1996; Gerin et al. 1997; Mathur & Nair 1997; Chengalur, 
de Bruyn & Narasimha 1999). In addition, a time-delay 
of 26±5 days between the two images of the lensed core 
has been measured (Lovell et al. 1998; see also Wiklind 
& AUoin 2002). New monitoring data in hand is expected 
to yield an improved measurement of this delay (J. Winn, 
private communications). The structure of the Einstein 
ring can be used to place constraints on the mass dis- 
tribution of the lens (Kochanek & Narayan 1992; Nair, 
Narasimha, & Rao 1993; Wikhnd & AUoin 2002). Were it 
not that PKS 1830-211 is seen through the Galactic disk 
(i.e. I — -f 12°, h — —5°) - leading to considerable optical 
extinction {Ai — 0.8 mag) and a field crowed by stars (e.g. 
Winn et al. 2002; Courbin et al. 2002) - this lens system 
would be well-suited to determine the Hubble Constant. 
Due to the severe extinction it has proven difficult to pre- 
cisely determine the center of the lens galaxy. A number 
of mutually inconsistent lens centers have been published 
over the years (e.g. Kochanek & Narayan 1992; Nair et al. 
1993; Lehar et al. 2000). Most recently, using the same 
HST observations, two different interpretations of the lens 
center were published (Winn et al. 2002; Courbin et al. 



2002). The inability to measure the lens center with the 
currently available observations severely limits the use of 
PKS 1830-211 to determine Hq. Forthcoming HST-ACS 
observations, however, will be able to settle this question. 

Besides being interesting because of Hq , the lens galaxy 
in PKS 1830-211 is the highest-redshift spiral lens galaxy 
known. The brightness of the lensed source has allowed 
a unique study of the molecular (e.g. CO and OH) and 
atomic (i.e. HI) interstellar medium (ISM) content (Wik- 
hnd & Combes 1996; Gerin et al. 1997; Mathur & Nair 
1997; Chengalur, de Bruyn & Narasimha 1999) in a spiral 
galaxy as seen KiS hj^ Gyr ago. Although only a single 
galaxy, it provides a test case for similar systems antici- 
pated to be discovered with planned instruments such as 
ALMA, LOFAR and SKA. 

The time base-line covered by these galaxies provide 
enough leverage to study the evolution of the mass- 
distribution and ISM of spiral galaxies in the redshift- 
regime where most of their evolution is expected to occur 
and also connects up to studies of damped Ly-a systems 
(e.g. Kauffmann 1996; Mo, Mao & White 1998). 

In this letter, we make a first attempt to model the kine- 
matics of HI in this high-redshift spiral galaxy, improving 
upon a previous more simple analysis by Chengalur et al. 
(1999). One of the initial aims, i.e. the determination of 
an independent kinematic lens center, has proven too dif- 
ficult given the current data, but constraints on the radial 
distribution of HI, its column density and kinematics are 
obtained. 
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used with 256 channels, giving 3.9 kms~^ per channel. 
The data were Hanning smoothed to give a resolution at 
halfwidth of 7.9 km s~^. The resulting normalized line pro- 
file is shown in Fig. 2. Errors are estimated from the rms 
scatter between channels at the baseline level well away 
from the line profile. 

3. KINEMATIC MODEL 

We model the HI absorption line, assuming a thin ax- 
isymmetric HI disk with an optical depth r(r) = tq x 
(r/ro)'^. The disk has a kinematic center {xk,yk) on the 
sky, an inclination i and position angle 9 (measured north 
to east on the sky) . The HI gas has a constant circular ve- 
locity Vc, as expected for an isothermal mass distribution 
of the lens galaxy and a systemic velocity Vs- The latter 
allows for a small velocity difference with respect to the 
redshift determined from molecular lines. For each sur- 
face brightness element T,^ (x) of the lensed radio source 
[x = {x,y) are coordinates on the sky], we then calcu- 
late (i) the radial component Vr{x) of the velocity vec- 
tor of the HI gas, where the line-of-sight pierces the disk, 
and (ii) the optical depth t{x) for the sight line normal 
to the disk, which is then multiplied by cos(i)~^ to in- 
clude the effect of inclination. We assume that the HI 
absorption line is isotropically (Gaussian) broadened by 
Chi (=FWHM/2.35), constant as a function of radius. 

We re-normalize the flux-density of the source, 
j T,v(x)(^x, to unity. The observed line profile becomes 
a convolution of a Gaussian velocity profile with the un- 
broadened line-profile^, i.e.: 



I{v) = {< 




S„(f) [l-r(f)] dx) (1) 



1.— The MERLIN 1.4-GHz continuum maps of PKS 1830- 
211 (North is up and East is left). Contours start at 25 mjy/beam 
and increase by factors of two. Overplotted is the best-fit kinematic 
model for the lens center of Nair et al. (1993). The two crosses 
indate the lensed cores A (NE) and B (SW). The thin elliptical 
contours indicate constant optical depth increasing by factors of 2 
outward. The optical depth contour passing through image B has 
r=0.064. 

2. MERLIN & WSRT OBSERVATIONS 

MERLIN long-track observations of PKS1830-211 were 
obtained at 1.4, 1.7 and 4.9 GHz on, respectively, March 
26, 25 and 8 in 1999. The data were flux and phase cal- 
ibrated at the telescope. We subsequently map the data 
using DIFMAP (Shepherd 1997) following standard pro- 
cedures, but perform additional self-calibration in the case 
that small phase and flux calibration errors are still present 
in the data. The resulting map at 1.4 GHz is shown in 
Fig.l. All maps will be presented in more detail in forth- 
coming publications. The 1.4-GHz map is used as a rep- 
resentation of the surface brightness distribution against 
which we observe the HI distribution in the lens galaxy. 
We use it to calculate the HI absorption line. We note that 
absorption really occurs at 1.42/(1 -I- zi) = 0.75 GHz, for 
which we do not have a radio map. However, the difference 
between the 1.4 and 1.7 GHz maps are small and we ex- 
pect the radio map at 0.75 GHz to be quite similar to that 
at 1.4 GHz. Using the 1.7 and 1.4 GHz radio maps to con- 
struct an estimated 0.75-GHz radio map is difficult given 
the different UV sampling. We therefore opt to solely use 
the 1.4-GHz map. 

The multi-frequency frontends of the WSRT were used 
to obtain a high signal-to-noise profile of the HI absorp- 
tion line at a frequency of 753 MHz. Compared to the 
previous observations of Chengalur et al. (1999) the data 
have superior spectral resolution and wider bandwidth as 
well as a better spectral baseline. In total we obtained 
5.5 hrs of data on PKS1830-211 with 12 MFFE's on Oc- 
tober 25, 1998. A backend bandwidth of 2.5 MHZ was 

^ This is correct only if the cthi is not a function of x and note also that the Gaussian is not normalized. 

^ Clean components are delta functions with finite flux that, convolved with a Gaussian beam, reproduce the radio map as presented in Fig.l. 
In total our model of the source brightness distribution includes 350 clean components. 

The line-width cthi might increase slightly because of unaccounted for minor velocity gradients between clean components. 



which is unity in the case of no absorption. Images of the 
radio sources - even though they have a continuous map- 
ping and surface brightness distribution - are always af- 
fected by the finite beam size. Hence, instead of using the 
MERLIN image of PKS1830-211 to model E„(x), we in- 
stead use the clean components^ . In this case the integral 
in Eq. 1 is simply replaced by a sum over the clean compo- 
nents. The resulting absorption line profile should there- 
fore more closely resemble the underlying line profile'^ for 
the deconvolved radio image of PKS1830— 211. The free 
parameters in our kinematic model are {i, 9, Vg, chi, tq, 7}. 
We furthermore set Vc = \/2 x ctdm — 266 km s^^, which 
follows from the isothermal lens mass model (i.e. a dark- 
matter dispersion of (TDM = 188kms^^; in agreement with 
Winn et al. 2002), ro = ^.'5 (arbitrary value but roughly 
the Einstein radius), and the kinematic center is set to the 
lens centers inferred from the different authors discussed 
in Sect. 1. We register the optical and radio data, as- 
sociating the bright radio cores with the optical quasar 
images. 

4. RESULTS 

Constraints on the kinematic model are the follow- 
ing: (i) the WSRT HI hne-profile with 238 channels with 
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3.9 kms~^ width each, with an average l-cr measurement 
error of 0.00147 (Fig. 2; three channels are removed due 
to band-pass edge uncertainties), and (ii) the velocity dif- 
ference between components A and B of 147(±l)kms~^ 
(Wiklind & Combes 1998). We minimize ^ a-s com- 
monly defined - using a Downhill Simplex method with 
simulated annealing (Press et al. 1992), varying all six 
free parameters simultaneously. 

We only explore models where the kinematic centers 
have been fixed at positions (notably often mutually ex- 
clusive) as previously published in the literature. Given 
degeneracies in the kinematic model, the current spatially 
unresolved line-profile and potential other issues (see be- 
low) , we feel that a more sophisticated kinematic model is 
at this point not warranted. 

The results are listed in Table 1. A number of conclu- 
sions can be drawn: 

(1) The inclination of the HI disk, i = {17 - 32)° ± 2° 

(1-cr), is robust and nearly independent from the assumed 
kinematic center of the lens galaxy. This confirms that the 
lens galaxy is most likely a close-to face-on spiral galaxy, 
as previously proposed (e.g. Winn et al. 2002). 

(2) The position angle of the disk varies between = -|-34° 
and —15°. This result is clearly less robust than the incli- 
nation angle, due to degeneracies with the inclination and 
kinematic center. 

(3) The line-width, am = (39 - 48) ± 1 kms'S is also 
robust. The value appears high compared with local spi- 
ral galaxies and might be the result of increased star for- 
mation, stirring up the HI gas to higher turbulent veloci- 
ties. The line-width of molecular gas is also relatively high, 
~30 km s^^ (Wiklind & Combes 1996). Even though lower 
than what we infer, such difference can be expected be- 
tween colder molecular gas, presumably in clouds, and HI 
in the atomic (versus molecular) phase. However, the use 
of clean components to represent the, in reality, smooth 
Einstein Ring could somewhat broaden the fitted line- 
width because they undersample the velocity field (i.e. the 
field in between clean components will shows a velocity 
gradient that is not sampled). The spectral resolution is 
negligible. 

(4) The optical depth varies considerably as a function of 
radius for some of the models, with typical l-cr statisti- 
cal errors on tq of ±0.01. In particular the models with 
the lens centers from Winn et al. (2002), Kochanek & 
Narayan (1992) and Courbin et al. (2002; only for model 
Sp) seem to prefer steeply rising optical depths as a func- 
tion of radius (i.e. 7^0). The other models also prefer 
rising optical depths, but more moderately. All models 
give a considerable HI depletion toward the lens center, 
as is also seen in some local galaxies, which could be the 
result of strong past and ongoing star formation (see also 
-3- above). This depletion can already been inferred from 
the fact that, even though image B is closer the inferred 
lens centers, it experiences less HI absorption (Fig. 2). 

(5) Plotting r(r), we find that all models give r « 0.085 
at r = 5.0 h^Q kpc along the disk, except for the lens 
center (Sp) from Courbin et al. (2002) for which 0.06 
is found. However, from a lensing perspective it is un- 
likely that the latter position is that of the dominant 



lens galaxy. The HI column density is A'hi = 4.6 x 
102i(rs/100 K)(crHi/kms-i)T cm'^ (given our defini- 
tion of the line profile in Eq. 1). Assuming a "canonical" 
value of the unknown spin temperature of Tg = 100 K, 
as in Chengalur et al. (1999), r w 0.085 at a radius 
of 5hjQ kpc in the disk and am « 44 kms~^, we find 
A^Hi ~ 2 X 10^^ cm~^, similar to what was already found 
by Chengalur et al. and also similar to that of local spiral 
galaxies and damped Lyo: systems. 

1 1 1 1 1 1 .1 1 1 1 1 1 1 1 



>> 


1.000 'i 


DensI 


0.990 H 


Flux 


0.980 H 


sed 


0.970 'r 


Normali 


0.960 'r 

0.950 r 



■ I ■ I ■ ' 

-300 -200 -100 100 200 
Velocity 

Fig. 2.— The WSRT Hl-linc profile ovcrplotted with the best-fit 
kinematic model from Fig.l. Errors are 0.00147 in normalised units. 
The velocity is in units of kms~^. A and B indicate the velocity 
positions of the two lensed cores of source. 

We note that the range of parameter values are typically 
larger than their 1 a errors, since the former reflects the 
change as function of chosen kinematic center. 

The results have to be seen in light of some additional 
cautionary notes: First, we see that the values are 
~2 times the DOF=229. This is mostly due to a mismatch 
between the model and the observed HI line around A and 
B (see Fig. 2). The lensed core of the source, however, is 
known to be time- variable (Lovell et al. 1998). Since, the 
MERLIN and WSRT data set were not obtained at ex- 
actly the same epoch (5 months in between), this might 
explain the mismatch. WSRT data from difi^erent epochs 
(not presented here) indeed show some minor differences 
between the observed HI lines. 

Second, our models are azimuthally symmetric and as- 
sume an HI covering factor of 100%. If the HI column 
density (i.e. optical depth) is a function of angle along the 
disk, it could change the inferred kinematic results. Even 
though we believe the effect to be small, it can not fully be 
excluded. A spatially resolved HI absorption line profile 
would help to eliminate this possibility. 

Third, we have assumed the 1.4-GHz maps to be a rea- 
sonable representation of the 0.75-GHz continuum struc- 
ture of the lensed source, as discussed previously. If the 
spectral index is a strong function along the lensed source, 
it might also affect the results. 

Fourth, we have assumed that the rotation curve of the 
spiral galaxy is flat. This is consistent with observations 
of local spiral galaxies with similar rotation velocities, but 
might not be a correct assumption for a galaxy at 2: = 0.89. 
Again, a spatially resolved HI absorption line profile would 
be able to test this assumption. 
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1.73 


0.049 


449 


Courbin et al. (2002), Sp 



Table 1 — Best-fits parameters of the kinematic models (see Sect. 3-4, also for errors). 



5. DISCUSSION & CONCLUSIONS 

We have modeled the neutral hydrogen (HI) distribu- 
tion and kinematics in the spiral galaxy of the lens sys- 
tem PKS 1830-211 at z=0.89, using a MERLIN 1.4-GHz 
continuum map of the Einstein ring and an WSRT HI ab- 
sorption line. The main numerical results are summarized 
in Sect. 4. 

In general, we find that the HI disk is seen closc-to face- 
on, as previously suggested (e.g. Winn ct al. 2002). We 
find evidence for HI depletion toward the galaxy center and 
a column density of TVhi w 2 x 10^"*^ cm~^ at 5 h^^ kpc, 
comparable to those of low redshift spiral galaxies. Our 
data are insufficient to discriminate between the different 
lens centers as published by different authors and we have 
to await new HST-ACS observations. 

We emphasize some cautionary notes, listed in Sect. 4, 
regarding our results. Even so, we believe that our anal- 
ysis is a good first attempt to constrain the kinematics of 
neutral hydrogen in a spiral galaxy at relatively high red- 
shift, within the current observational limits, and should 
be regarded as such. 

Our analysis shows that the combination of HI observa- 
tions with gravitational lensing - in particular radio Ein- 
stein rings - potentially provide a strong tool to study 
spiral galaxies to high redshifts. A similar program to 
combine lensing and stellar kinematics of E/SO galaxies 
out to z ~ 1 is currently underway (i.e. the Lenses Struc- 
ture & Dynamics (LSD) Survey; e.g. Koopmans & Treu 
2002, 2003; Treu & Koopmans 2002). 

We anticipate that spatially resolved HI absorption line 



profiles and radio continuum images at the rcdshiftcd HI- 
line frequency can be obtained of many similar systems 
in the future with e.g. ALMA, LOFAR and SKA. This 
would resolve many of the issues outlined in the paper and 
provide us with an unprecedented probe into the internal 
structure of spiral galaxies at high redshifts. By probing 
spiral lenses over a range of redshifts, their evolution can 
be studied (e.g. Koopmans & de Bruyn 1999), provid- 
ing a complementary tool to study spiral galaxies between 
those at low redshifts using HI emission and damped Ly-a 
systems at very high redshifts. 

We note, however, that none of the other studies can 
provide two independent constraints, i.e. from lensing and 
kinematics, on the same mass distribution of the spiral 
galaxy. Their combination can break degeneracies inher- 
ent to each technique separately. Our results, admittedly 
based on a single system, show that very interesting con- 
straints can be obtained for spiral galaxies at high red- 
shifts. 
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